Statistics of Spectra for One-dimensional Quasi-Periodic Systems at the 

Metal-Insulator Transition 
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We study spectral statistics of one-dimensional quasi-periodic systems at the metal-insulator 
transition. Several types of spectral statistics are observed at the critical points, lines, and region. 
On the critical lines, we find the bandwidth distribution Pb (id) around the origin (in the tail) to have 
the form of Pb(vj) ~ w a (P B (w) ~ e' 13 ^) (a,/3,7 > 0), while in the critical region Pb(w) ~ w a 
(a' > 0). We also find the level spacing distribution to follow an inverse power law Pg(s) ~ s~ s 
(5>0). 
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The universality of statistical properties of energy 
spectrum is ubiquitously found in quantum physics 
ranging from quantum chaos to quantum chromo- 
dynamics |3(- In condensed matter physics, the metal- 
insulator (MI) transition in disordered electron systems 
gives a notable example of such universality. In the 
metallic side, the energy level statistics can be described 
by random matrix theory Q . The level spacing distribu- 
tion is close to the Wigner surmise. The universality class 
is classified by the symmetry of the system. In the insu- 
lating side the level spacing distribution is Poissonian. At 
the critical point, new statistical features emerge 0.0 
They are different from random matrix theory and Pois- 
son statistics. The universality of these new statistics 
dubbed as "critical level statistics" is also classified by 
symmetry of the ensemble 0, H, 0, El • 

Quasi-periodic systems are other interesting ones 
showing the MI transition. The universality of the tran- 
sition has been characterized by multifractal structures 
of bandwidths and wavefunctions (see Ref. 0] for a re- 
view). The Harper model 0,0] is obtained by a gauge 
fixing of two-dimensional electrons on the square lattice 
in a uniform magnetic field. Especially incommensurate 
limits of the flux per plaquette, such as the inverse of the 
golden mean a = , have been extensively studied. 

In the atomic scale, the penetration of flux for each pla- 
quette requires an enormous strength of magnetic field, 
but nowadays it is realized in laboratory on a quantum 
dot lattice 0. 

The quasi-periodic system can be seen as a system be- 
tween "periodic" and "random," thus giving an interest- 
ing example of quantum "chaos." Classically an open 
orbit (separatrix) of equi-energetic curve appears at the 
MI transition point. But the quantum transition cannot 
be understood by the classical consideration. It is not 
characterized by a single energy band corresponding to 
the separatrix but the wavefunctions for all the energy 
bands become critical with power-law decay. The energy 



level statistics are different from disordered systems. It is 
not close to the Wigner surmise on the metallic side nor 
is the Poisson statistics seen on the insulating side (lif . 
Spectra of quasi-periodic systems have a fractal struc- 
ture 0], which hinders established techniques used in 
disordered case such as "unfolding" to extract universal 
properties of fluctuations. 

Recently Evangelou and Pichard investigated the 
bandwidth distribution Pb{w) (w : bandwidth) of the 
one-dimensional Harper model at the incommensurate 
limit of a. They found that it agrees with the semi- 
Poisson statistics Pb(w) — 4we~ 2w with the sub-Poisson 
number variance ^(E) ~ const. + \E (x < 1) at the 
critical point. They argued that it is related to the uni- 
versal characterization of quantum chaos. 

Although this suggests that a characterization by level- 
statistical idea is possible for the MI transitions of a 
larger class of quasi-periodic systems, the universal ap- 
pearance of the semi-Poisson statistics may be questioned 
from the experience in critical level statistics of disor- 
dered systems. In Ref.Q, it has been suggested the level 
spacing distribution may be described by a more gener- 
alized form P(s) = As a e~^ s (a, 0, 7 > 0), and it has 
been observed in numerical studies [lfj. A similar rich 
structure of critical level statistics should be expected for 
quasi-periodic systems. 

Motivated by this idea, we investigate a generalization 
of the Harper model, namely the extended Harper model, 
which is obtained from two dimensional electrons on the 
square lattice with next-nearest-neighbor hopping in a 
uniform magnetic field 0, 0] . We find that such a rich 
structure indeed emerges for the critical level statistics of 
this extended model, but with some substantial modifi- 
cations, even for the case treated in [l^. We explain this 
variety of critical level statistics by the quantum nature 
of the system. 

Let us first describe the model. We consider the 
two-dimensional fermions on the square lattice in a 
flux per plaquette tp with next-nearest-neighbor coupling 
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FIG. 1: (a) Transfer integrals of the extended Harper model 
(Eq.Q) and (b) its phase diagram (See H3). In region I the 
wavefunction (spectrum) is extended (absolutely continuous) 
and in region II it is localized (pure points). In region III and 
on the three boundary lines, it is critical (singular continu- 
ous). 
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FIG. 2: The bandwidth distribution at (A, y) = (2.0,0.4), 
and the one near the origin (inset). 



(Fig^a)). The Hamiltonian is given by 

H = cosp x + Xcosx + /icos(p x + x) + [icos(p x — x), (1) 

where p x — —lixify-^. The positive parameters, A and fi, 
are the strength of the kinetic terms for the perpendicular 
and next-nearest-neighbor directions respectively, while 
the horizontal one is normalized to unity. We take the 
Landau gauge ^{x n ,y n ) = ip x (x n )e %vyn and obtain the 
Schrodinger equation for ip n = il) x {x n ) as 



1 + // en-, j 2n(n + 



4>n 



+ i 



-A cos(2Tmip + v)ip v 



I + // cos j 2n{n —-)</?- 



4> n -i = Eif> n . (2) 



At fi = 0, this is reduced to the Harper model. For 
a rational tp = p/q where p and q are relatively prime, 
Eq. (J2J reduces to a q x q matrix with the boundary con- 
dition ipn+q — exp(ikq)tp n , which can be transformed 
to be tridiagonal 18] . To investigate the irrational 
limit of ip — > a, we use the rational approximation of 
ifj = Fj-i/Fj where Fj is the Fibonacci number satisfy- 
ing F j+1 = Fj + Fj_ x with F = 1 and F x = 1. 

The phase diagram is shown in Fig^b). The scaling 
property of the bandwidth was studied in 01 ■ For large 
q, the minimum and maximum are obtained at either k = 
or ir/q for fixed v ^^|. We numerically diagonalized the 
matrices there setting v = for simplicity, and obtained 
the bandwidth distribution Pb{w). The normalizations 
are J Pb{w)cIw — 1 and (w) — L wPb(w)cIw = 1. 

First, we study the bandwidth distribution along the 
critical line A = 2. In Fig[3 we plot the Pb{w) at 
(A, n) = (2.0, 0.4) for j = 25, 27, 28. It shows a good 
convergence to a limit distribution, indicating the exis- 
tence of the limit of the bandwidth distribution at the 
incommensurate flux tp. For < /i < 1 , we find Pg(w) 
exhibits a continuous change. It has been concluded |l7| 
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FIG. 3: /i-dependence of index j3 along the critical line A = 2. 



that the scaling property of the bandwidths is invariant 
on this line, implying the system belongs to the same uni- 
versality class. The continuous change means that Pg (w) 
is not a universal quantity. 

In Fig |2 one sees linear behaviors at large w, implying 
an asymptotic form 



Pb(w) 



-/3w 



as w 



(3) 



where (3 > 0. The optimized values of /3 are shown in 
Table Q] It suggests a quadratic dependence of (3 on fi. 
The best fit is attained by a quadratic function ~ 1.4 — 
0.53/x 2 (Fig[3J). This rather simple dependence may be 
related to the universality class on this line. 

Pb(w) is estimated to be zero at the origin. Inset of 
Fig|21 shows Pb (w) near the origin, indicating a power 
law. To characterize this behavior, we make an ansatz of 



Pb(w) 







(4) 



where a > 0. The optimized values of a(/i) are shown in 
Table HJ One sees that ev(/x)s are stable, suggesting that 
a(/i) is independent of /j, at the incommensurate limit. 

These numerical results suggest that the overall distri- 
bution is given by a generalized semi-Poisson form 



P B {w) = Aw a e 



-f3w 



(5) 



3 



bo 
o 



1 

0.5 

-0.5 

-1 
-1.5 

-2 
-2.5 

-3 
-3.5 

-4 



oo 

'O 

b -q 



o 



3 
2 

£ ° 

g>-2 
- -3 

3 -4 

Qq -5 

00. 




o 



2 -1.5 -1 -0.5 0.5 1 



00-. 



log w 



Oo»Q.. 



q=17711 O 

P B (w) ~ exp(-1.93w 064 ) 



OOOqo 



00^ 



"doo°o - 



0.5 1 1.5 2 2.5 3 3.5 4 
log w 

FIG. 4: The bandwidth distribution at (A,/i) = (1.0,1.0), 
and the one near the origin (inset). 



where A is a normalization constant. The exact semi- 
Poisson form i.e. a = 1, /3 = 2 has been reported at 
(A,/z) = (2.0,0.0) [16(. Actually the normalization con- 
ditions allow only one parameter a with f3 = a + 1 and 

A = .,, . It is clear that such relations are not con- 

r(a+i) 

sistcnt with our results as shown in Table |I| Especially 
the semi-Poisson statistics with a = 1 do not reproduce 
the index of the power-law behavior around the origin at 
fj, = 0. Using one-parameter fit by (JjJJ, we get a ~ 0.7 
to the overall distribution, which is much smaller than 
that estimated by (0J. For // ^ 0, we get similar de- 
viations. Thus we may conclude that the semi-Poisson 
form for Pb(w) is only an approximation to the overall 
distribution. 

Next, we investigate Pb(w) on the other critical lines. 
It turns out that the behaviors are different from those 
on the critical line A = 2. Fig 0] shows the Pg(io) at 
(A,//) = (1.0,1.0). The large w behaviors do not show 
an exponential decay, but a milder one. Thus we make a 
generalized ansatz of 



Pb{w) 



as w 



(6) 



with < 7 < 1. The optimized curve agrees well with 
the obtained data (FigQJ. The behavior near the origin 
is also found to follow a power law (inset of FigQJ. On 
the line fj, = 1, the index a seems to be constant around 
1.5. We find that the numerical data on the critical line 
A = 2/i can be also described by each ansatz Q and JBJ- 
The index a is stable around 1.5 as in Table ITT1 which is 
near the values on the critical line fj, — 1. 

The overall distribution on the critical lines [i — 1 and 
A = 2/i except for the bicritical point may be cast in a 
form 



This form of distribution was suggested in Ref.[7j for the 
level spacing distribution of disordered systems at the 
MI transition. The 7 ^ 1 case has been encountered in 
a 3D disordered system with orthogonal symmetry |l0| . 

/r(a±av" 

The normalization conditions constrain p 



and A = ^C|) 



P B {w) =Aw a e 



(7) 



Our analysis does not support these 

constraints. Again, we tried a one-parameter fit varying 
0(7) in (J7J with the observed value of 7(a) fixed. We got 
values much smaller (larger) than the observed indices a 
(7) around the origin (in the tail) for all the points we 
studied on the critical lines /i = 1 and A = 2/i. 

We also investigate the bicritical point at (A,/i) = 
(2.0, 1.0). However, the convergence of the obtained dis- 
tribution is slow, and we have not got information of the 
incommensurate limit. 

The values of a and 7 on the critical lines /i = 1 and 
A = 2/i are considerably smaller from those on A = 2 
(7=1 there). This smallness of a and 7 indicates the 
tendency of Pb(w) on these critical lines to broaden both 
to the origin and the tail compared to Pb(w) on A = 2. 
This may imply that the transport property becomes dif- 
ferent when /i gets sufficiently large. To supplement this 
observation, we investigate the bandwidth distribution 
in the critical region (the region III in Fig^Jb)). Pb{w) 
follows an inverse power law Pb(w) ~ w~ a (a' > 0) 
there (FiglSJ. The values of a' are summarized in Table 
[0 This behavior sharply contrasts with the one seen on 
the critical lines. The divergence of Pb(w) at the origin 
implies the dominance of relatively flat bands, indicating 
that in most of the bands the wavefunction is close to be 
localized. Also the power law decay in the tail means the 
appearance of bands whose wavefunction is close to be 
extended. 

From the two-dimensional point of view, the next- 
nearest-neighbor hopping term (/i dependent term) is 
dominant in Eq.Q in the critical region. Since y- 
translation is canonically conjugate to ^-translation in 
the uniform magnetic field, the extension (localization) 
of the wavefunction in //-direction must be balanced by 
the localization (extension) in x-direction to satisfy the 
uncertainty principle. When [i is small, ^-dependence of 
the wavefunction is determined by A, the anisotropic pa- 
rameter for //-direction. On the other hand, when fi is 
sufficiently large, as fj, acts on x- and //-directions equally 
the wavefunction is sensitive to fj, by both of its x- and 
//-dependences. This makes the MI transition in the crit- 
ical region different from that on the critical line A = 2, 
resulting in the increase of the bands close to be local- 
ized as well as the bands close to be extended. Thus the 
appearance of the inverse power law in the critical region 
and the differences of a and 7 among the critical lines 
are consequences of quantum nature of the system. 

We also investigate the gap distribution Pg(s)- The 
gap distribution at (A, /1) = (2.0, 0.0) has been known to 
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FIG. 5: The band width distribution and level spacing dis- 
tribution (inset) at (A, /i) = (2.0,2.0). 



follow an inverse power law |l9l l'2(| , which diverges at 
the origin 

-S 



Pg(s) 



(8) 



with 5 ~ 1.5. We find Pg(s) along the three critical lines 
and in the critical region also follows an inverse power 
law (inset of Fig [5}. The values of S are shown in Table 
□ and OH 

It would be interesting to examine a similar statistical 
distribution for the wavefunctions. The characterization 
through level-statistical idea may have a good possibility. 

In this letter, we have seen that a systematic charac- 
terization by a level-statistical idea for the quasi-periodic 
system is possible. We have investigated the bandwidth 
Pb(w) and gap distributions Pg( s ) for one-dimensional 
quasi-periodic Schrodinger equations at the MI transi- 
tion and confirmed their variety as expected. We found 
a power law Pb(w) ~ w a near the origin and a general- 



ized exponential decay Pb(w) 



-pw~> 



(a, (3, 7 > 0) on 



the critical lines, while in the critical region Pb(w) fol- 
lows an inverse power law ~ w~ a (a' > 0). We gave the 
explanation for this variety of Pb (w) by quantum nature 
of the system. A surmise of a form Pb(w) = Aw a e~^ w ' 1 
gives only an approximation. The gap distribution Pg(s) 



shows an inverse power law 



,-5 



for whole the phase 



diagram. For the bicritical point, we did not get a con- 
clusive result. 
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